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"st" 1 Abstract 

This talk summarizes recent progress in lattice QCD for dense quark matter. 
' The emphasis is on the insights obtained from analytical results derived within 

■ chiral perturbation theory. 

> '■ 

^ 1 Introduction 

Despite the fact that all observed processes conserve baryon number, the total 
number of baryons present today is far greater than the number of anti-baryons. 
It is a long standing challenge to explain this excess of baryons. If we accept 
the baryon imbalance as a fact and ask: 

What is the preferred ground state of strongly interacting matter given a 
specific baryon density and temperature ? 

we are faced with an equally long standing challenge: The phases of strongly 
interacting matter must be determined by non-perturbative means, and the 
only known first principle method, lattice QCD, is extremely challenging to 
apply when there is an imbalance between baryons and anti-baryons. 
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Even though these two challenges are closely related their solutions will 
presumable be of a very different nature. To understand the generation of the 
baryon number most likely requires physics beyond the standard model. To 
understand the implications of a non zero baryon number on the ground state 
of strongly interacting matter requires the development of new non-perturbative 
tools for QCD. 

This talk focuses on the latter challenge so let us explain in detail what 
we are up against when we seek to describe baryon asymmetric matter from 
first principles (even when we are modest, in that we do not aim to explain 
where this asymmetry originally came from, but simply want to adopt this into 
QCD). For simplicity we will also ignore electromagnetism. 

Rather than fixing the density the natural approach is to start with the 
Grand Canonical partition function, Z, where we fix the quark chemical poten- 
tial, and then determine the quark number n (the baryon number is 3n). 
The average quark number is 

(n)=^logZ(/i). (1) 

It is straightforward to include \x in the partition function. Since the chemical 
potential is conjugate to the quark number, the combination that appears in 
the Lagrangian is 

/in = [iq^q = nq^q. (2) 
The simple extension of the Dirac operator in the Lagrangian of QCD 

£qcd = q{Dnlr, + M7o + m)q + Gluons (3) 

thus accommodates the chemical potential into QCD. 

It is well known how to implement this extension on the lattice [H [2] . Since 
\jl enters as the zeroth component of a vector potential the choice with the 
correct continuum limit is [T] 

^Lattice QCD = - + e afl q xl0 U xx+(> q x+6 + e~ a % + o7o^ ia:+5 <fa + ••• (4) 

For two color QCD, when the fermion determinant is real, lattice simulations 
[31 HI El [7] based on this discretization give a rich phase diagram that agrees 
with theoretical expectations [HEIITU]. For three colors, however, the fermion 
determinant at non zero chemical potential is not real and positive 

det 2 {D + ^70 + m) = \det(D + u 7o + m)\ 2 e 2ie . (5) 

While there is nothing wrong physically with a complex valued fermion deter- 
minant it is a major problem for lattice QCD. The Monte Carlo sampling of 
gauge field configurations, which is the back bone of lattice QCD, only works 
if the measure in the Euclidean partition function is real and positive. This is 
the QCD sign problem: at fi ^ direct Monte Carlo sampling is not possible 
and standard lattice QCD breaks down. 
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Several methods have been engineered to circumvent the sign problem, see 
table [H Each have given important first principle insights into QCD at non 
zero and a consistent picture emerges at small fJ,/T. All methods face serious 
challenges more or less directly induced by the sign problem. In order to deal 
with these challenges analytical results are of utmost importance. For example, 
in the method of Ejiri [28] the distribution of the phase of the fermion determi- 
nant is assumed to be Gaussian. Here we show that the Gaussian form comes 
out analytically for small chemical potentials, thus confirming the assumptions 
of this lattice method. However, as we also show, the distribution changes to 
a Lorentzian form for larger jj,. This computation is carried out within chi- 
ral perturbation theory, which is the low energy effective theory for QCD in 
the chirally broken phase. The analytic results in this way both justifies the 
assumptions on which the method of Ejiri is based and show its limitations. 

Chiral perturbation theory also allows us to understand the distribution of 
the values which the Euclidean quark number operator assumes over the gauge 
fields. The results not only give us a direct insight in the way the total baryon 
number forms, it also shows how the Complex Langevin method can deal with 
the sign problem. 



2 Fixed phase of the determinant 

Perhaps the most direct way to understand the severity of the QCD sign prob- 
lem is to consider the distribution of the phase of the fermion determinant 

(5(9 - 9 )) 1+1 d8 /dA | det(2? + A17b + m) |2 e 2^ e -SvM d °- W 

The 5 function allows us to rewrite the unquenched distribution in terms of the 
phase quenched one times a phase factor and a normalization 

(5(9 - 0')>i+i = e 2ie ^(5(9 - 0')>i+i* • (7) 
^i+i 

The complex nature of the unquenched 9 distribution is typical of the sign 
problem. Here it takes the simplest possible form. Nevertheless, the effect 
of the phase is dramatic: it leads to exponentially large cancellations in the 
volume. This follows directly since Z\+\* / 'Zi+i ~ exp(V). These dramatic 
cancellations are what makes the sign problem severe. 

To understand how the cancellations take part we have computed the 9 
distributions to leading order in chiral perturbation theory. For /j, < m n /2 
with the phase angle 9 £ [— tt, tt] we have [38j 

-. oo 

(5(9-9')) 1+1 = 1 P *fl+AGo £ e -(s+a nff )VAG 0j (8) 
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Method 


Idea 


Challenge 


Reweighting 

rm [T21 mi 

|A±i LL^4 LiiiJ 


Absorb the sign in the observable 


Exponential cancellations 


Taylor expansion 
[T4~l fTHl [TBI n~7l fTRl 


Expand at [i = 


Higher order terms 


Imaginary /x 
PTOl [201 12TI \T2l 


Determine the analytic function 


Control the extrapolation 


Density of states 
1231 [231 [231 [261 1271 [281 


Use the distribution of the phase 


Determine the distribution 


Canonical ensemble 
[291 EH ED E21 E3] 


Work at fixed baryon number 


Fix the baryon number 


Complex Langevin 

[311 E3 ESI EH 


Stochastic flow in complex plane 


Make it work for QCD 



Table 1: Summary of the main methods used to circumvent the sign problem, the 
main idea used and the main challenge faced when using this method. 



where 



Vm 2 T 2 ^ Ko( m -" n ) 2im 
AG = l2^V^i^ (00B h(«-l). (9) 



n 2 T 



Notice that AGo is the difference of the free energy in the full and the phase 
quenched theory 

|l±L = e -A Go _ (1Q) 

Z1+1* 

The 9 distribution takes the form of a Gaussian modulo 2ir. This is in agreement 
with observations on the lattice [28] and is the natural form suggested by the 
central limit theorem. However, at larger values of the chemical potential the 
quark mass enters the support of the Dirac spectrum [391 SO] and fluctuations 
of the phase will be far greater [H I H2l 1351 133] . 

To leading order in chiral perturbation theory for [i > m n /2 the 9 distribu- 
tion is a Lorentzian modulo 2ir times the phase factor |43| : 



/xoa ottw ^ VLb smh(VL B ) 

(d{29 - 26 i+i = e -— -, (11) 

7r cosh(l/ Lb) — cos(2#) 
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Figure 1: The quenched distribution of the phase of the fermion determinant. The 
dashed curves are standard Gaussian (left) and Lorentzian (right). The full curves 
are the corresponding distributions modulo 2ir and tc relevant for jj, < m n /2 and 
fi > respectively. 



where 

z,,, 

(12) 



Zi+i VL 



Zi+i* 

This demonstrates that the assumptions of the central limit theorem are not 
satisfied and that the method of Ejiri must be examined carefully for fi > m^/2. 

The Gaussian form (for /i < m n /2) and the Lorentzian form (for [x < m^/2) 
modulo 2tt are illustrated in Fig. [Tj 

3 Distribution of the quark number over 6 

Since chiral perturbation theory deals with pions alone, the quark number 
vanishes (we do not consider the effect of pion fields with non-trivial topological 
winding). However, the distribution of the quark number over 9 is non trivial 
even when evaluated within chiral perturbation theory. The form it takes gives 
us direct insight into the noise produced by the pions in lattice QCD. 

With fx < m-n/2 we find to leading order in chiral perturbation theory |43j 

oo Pin 2i8 

(n8{e-e>)) 1+l = v I E (i + ^^f^) 7 ^ r ^ +2 ™ )2/AGo+AG ° ) (i3) 

n=— oo V ^ 

where 

I /j=flim-^AG (-^/i)V (14) 

This distribution gives the quark number measured on a set of configurations 
for which the phase of the determinant is between 9 and 9 + d9. 
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Note that the normalization is by the two flavor partition function, so that 
the total quark number is the integral over the phase angle. 

For n > m-Tr/2 the distribution again changes its form drastically [43j 

(n 6(26 - 29>)) 1+1 = -^^]e^e^ __L_ , (15) 

Prior to discussing the importance of these results let us present the quark 
number distribution as predicted by chiral perturbation theory. 



4 Distribution of the quark number 

The distribution of the values which the quark number operator, 

n(/x)^Tr 70 , (16) 
D + fijo + m 

assumes as the gauge fields vary is also accessible within chiral perturbation 
theory. 
Since 

n (uY = ( Tr ^ ) * = Tr ^ = -n(-u), (17) 

the quark number operator is in general complex (it is purely imaginary at 
H = 0). The fluctuations of the quark number thus occur in the complex plane. 
The distribution of the quark number in the complex plane is by definition 

y) = (s(x- i(n(/i) - n(-/x))J 8 (y + i-(n(p) + n(-^ \ . 

1+ \l8) 

In order to evaluate this average within chiral perturbation theory both 5 func- 
tions are expanded in terms of the moments of the real and imaginary parts of 
the quark number. The first moment, i.e. the average quark number, vanishes. 
The higher moments are off-diagonal susceptibilities and are non-trivial even 
within chiral perturbation theory |44j . 

For n < m-x/2 we find to 1-loop order [44] 



P 1+1 (x y) = - e -( x -»i) 2 /(xld+Xud) e^y+vi) 2 /(xid-Xud) (19) 

^ixLY - {x B ud ? 



where 



d 2 

= , , AG (fj,x, fi 2 ) 



(20) 

Ml=M2=M 



d 2 

x ud = , , AG (-pi,p 2 ) 
dfiidfi 2 
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(21) 

/il=M2=M 



Note that the distribution factorizes into the distribution of the real and imag- 
inary part of the quark number. 

For ji > m n /2 the higher moments diverge and lead to power law tails of 
the distribution of the quark number operator [44] . In the quenched case, for 
example, the distribution has a inverse cubic tail. 

In the confined phase the pions dominate the free energy and the quark 
number is thus generated by subleading terms. A central aspect of the sign 
problem is to pick up this small baryon signal from the background produced 
by the pions. The results presented above give the analytical form of the pion 
background and thus the foundation for extracting the quark number. 



5 Cancellations 

The non trivial results for the distribution of the quark number with 9 and 
the distribution of the quark number itself should of course be consistent with 
zero average quark number. Therefore, upon integration over 9 and over the 
complex quark number plane respectively, we must find a vanishing average 
quark number. Here we show analytically how these cancellations occur. 



5.1 Integration over 6 

The total quark number is obtained from its distribution over 9 by integration 
over the phase angle 

(n) 1+1 = / d9(n6(9-9')) l+1 . (22) 



We now discuss how this becomes zero within chiral perturbation theory. 

First we consider the case fi < m n /2. Here we can write the distribution of 
the quark number over 9 as a total deriative of the 9 distribution times vj 

oo 1 J lid 

<„ M > 1+1 - „ £ IA^e-<««»>7™ ( 23) 

n=— oo v u 

= vr~Me-e')) 1+l . (24) 

1% a9 

From this we see that all phase angles, 9, are essential in order to obtain the 
desired vanishing value. Moreover, we see that these exact cancellations damp 
the exponentially large amplitude of the distribution completely (recall that 
AGq is extensive). Therefore, unless the 9 integration is carried out exactly 
one will find an exponentially large contribution to the quark number from the 
pions. 

Also for [i > m-x/2 we can write the distribution of the baryon number 
over 9 as a total derivative. However, this time it is not the derivative of the 
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^-distribution but rather 



(n <5(20 - 200)1+1 = - 



7T L /J, J 



1 rFL/ 




log(e 




(25) 



Again this implies that all angles contribute in an essential way to the full 
quark number. The need for an exact integration over 9 is even more urgent for 
\x > m n /2, since the amplitude now grows exponentially fast with the volume 
even at mean field level. 

5.2 Integration over the complex n plane 

Given the distribution (|19p of the baryon number, P^ +1 (x,y), we obtain the 
average quark number upon integration over the complex quark number plane 



The analytical result for the distribution shows precisely how the vanishing 
expectation value of the quark number is realized 



We see that the total quark number is obtained only after a detailed cancella- 
tion between the contribution from the real part and the imaginary part. The 
challenge for unquenched lattice QCD is to identify the baryon signal within 
these massive cancellations due to the pions. The analytical result from chi- 
ral perturbation theory shows precisely how this may be realized within the 
framework of the Complex Langevin method, for a detailed discussion see |44j . 



The description of the macroscopic phases of strongly interacting matter based 
directly on the microscopic QCD dynamics is hampered by the QCD sign prob- 
lem. 

Despite the complex nature of the sign problem substantial progress has 
been made recently. We have derived the distribution of the complex phase of 
the Dirac operator analytically and the distribution of the quark number over 
the phase, as well as the distribution of the quark number itself. These results 
give direct analytical insights into the numerical challenges faced by the density 
of states method and the Complex Langevin method when applied to QCD at 
non zero chemical potential. 

The new results form a novel branch of applications of chiral perturbation 
theory. This branch has its root in the evaluation of partially quenched averages 
by the replica method [35]. 




(26) 




(27) 



6 Summary 



s 
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